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5 ■ Abstract 

^J . We examine the cosmology of the two recently proposed scenarios for a five dimen- 

1-5 I sional universe with localized gravity. We find that the scenario with a non-compact 
fifth dimension is potentially viable, while the scenario which might solve the hierarchy 



^H. problem predicts a contracting universe, leading to a variety of cosmological problems. 
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The main theme of the first 20 years of the hierarchy problem has been to modify parti- 
cle physics around the TeV scale. In the past two years it has become evident that another 
successful route to solving the hierarchy problem is to modify the nature of gravitational 
interactions at distances shorter than a millimeter |jl], ^, |^. This modification can be most 
simply achieved by introducing compact extra dimensions. One recent proposal which at- 
tracted enormous attention is to lower the fundamental Planck scale M^ all the way to the 
TeV scale, by introducing large extra dimensions [^. The observed Planck scale is then just 
an effective scale valid for energies below the mass scale of the Kaluza-Klein (KK) excita- 
tions. The consequence of this proposal is that the necessary size of the extra dimensions is 
determined by the KK reduction formula 

where M* is the fundamental Planck scale of the order of 1 TeV, Mpi = 10^^ GeV, and n is 
the number of extra dimensions. Applying this formula for one extra dimension one obtains 
R ~ 10^^ cm, which would immediately suggest that this possibility is excluded because 
gravity would be modified at the scale of our solar system. For n > 2, R is sufficiently 
small so that this scenario is not excluded by short-distance gravitational measurements. 
However, Randall and Sundrum recently realized [^ Q that the case of one extra dimension 
is very special, and the naive KK reduction formula (|l]) may not be applicable in this case. 
The basic reason behind this is that in this scenario the standard model fields have to be 
confined to a three-dimensional wall ( "three-brane" ) , and such branes act like sources for 
gravity in the extra dimensions. The behavior of Green's functions in one dimension is 
dramatically different from the case of two or more dimensions. Indeed, in one dimension 
Green's functions grow linearly, while in the case of more than two dimensions there is an 
inverse power law l/r*^"^ (and logarithmic growth for d = 2). Thus branes act like small 
perturbations on the system for the case of two or more extra dimensions, and one expects 
the KK reduction formula (0) to be applicable. For just one extra dimension, however, the 
presence of branes can significantly alter the bulk gravitational fields, which may invalidate 
arguments based on the naive KK reduction formula. 

Indeed, Randall and Sundrum (RS) have presented a new static classical solution to 
Einstein's equations with one extra dimension (taken to be S^/Z2), and branes with non- 
vanishing tensions placed at the orbifold fixed points [Q. In this solution, for large brane 
separations, the effective four-dimensional Planck scale is independent of the size of the 
extra dimension, in agreement with the expectation that (|l]) is invalid for the case of one 
extra dimension. For their solution to work RS found that the brane tensions must be of 
opposite sign, and that there must be a negative bulk cosmological constant which stabilizes 
the system. 

Two possible applications of this solution have been proposed 0, 3 . In one case (which 
we refer to as RSI), our universe ("the visible brane") is the brane with negative tension, and 
the exponential "warp-factor" appearing in the RS solution will yield a natural new solution 
to the hierarchy problem ^. In the second (RS2) proposal p, the visible brane is the one 
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with positive tension. In this case the hierarchy problem is not solved, however, the second 
brane can be moved to infinity, thus providing an exciting example of a non-compact extra 
direction 0, which nevertheless correctly reproduces Newton's law on the visible brane. 

Binetruy, Deffayet and Langlois (BDL), however, have pointed out recently [0 that five 
dimensional theories with branes tend to have non-conventional cosmological solutions, once 
matter on the walls is included.]^ In this letter, we apply the results of BDL to the models 
presented by RS. We find that the equation for the scale factor on the visible brane (for small 
matter densities) coincides with the conventional Friedmann equation, up to the overall sign 
of the source terms. This sign depends on the sign of the cosmological constant (tension) on 
the visible brane. In the case of negative brane tension the source terms of the Friedmann- 
like equation have the opposite sign from standard cosmology. The change in sign implies 
that the universe would collapse on a timescale on the order of the Hubble time at the start 
of the expansion, for any matter with energy density p, pressure p, on our wall with an 
equation of state p = wp and w < 1/3. So whereas in the radiation-dominated phase the 
universe expands as in the conventional cosmology, after the transition from the radiation- 
dominated to matter-dominated (or quintessence-dominated) epoch when the universe was 
a few thousand years old, the universe would not expand as in conventional cosmology, but 
would rather collapse within a few thousand years. As we argue, this conclusion relies on 
knowing the expansion rate during the era of matter-radiation equality, which is provided 
by the success of the standard big-bang nucleosynthesis (BBN), and the current baryon and 
radiation densities. So in order to avoid this conclusion this scenario requires either a non- 
standard BBN or a modification of the RSI solution, for example through the introduction 
of additional fields. We emphasize that crucial to these conclusions is the fact that we live 
on a negative-tension brane, and that there is only one extra dimension. Both of these facts, 
however, are also crucial ingredients to the solution to the hierarchy problem presented by 
RS. For the case of the positive brane tension, however, the conventional expanding solution 
is reproduced in this model. 

We begin by summarizing the work of BDL [j^. The scenario considered here is a five- 
dimensional spacetime compactified on the line segment S^/Z2. The bulk coordinate is 
labeled by y, which is taken to be in the interval —1/2 < y < 1/2, where the points —1/2 
and 1/2 are identified. In this notation the coordinate y is dimensionless. The Z2 symmetry 
identifies the points y and —y, so we can restrict ourselves to < y < 1/2. Two three- 
branes are placed at the fixed points of the discrete symmetry, y = and y = 1/2. In our 
notation, the visible brane is located at y = 0, and a hidden brane is located at y = 1/2. The 
compactness of the extra dimension requires the existence of two branes, since each brane 
must absorb the gravitational fiux lines from the other. 

The most general metric for a five-dimensional spacetime which preserves three-dimen- 
sional rotational and translational invariance is given by 

ds^ = n^{T, y)dT^ - a^{T, y)dx^ - b^{T, y)dy^. (2) 

Note that here we use the metric signature (-|-,—,— ,—,—). The induced metric on the visible 
*For other recent results on the cosmological aspects of theories with large extra dimensions see [8-16]. 



brane is obtained by evaluating the metric tensor a.t y = 0. The metric is determined by 
solving Einstein's equations in the presence of some energy density in the bulk and on the 
three-branes. A solution which describes a static four- dimensional Lorentz invariant universe 
is given by a{T,y) = n{T,y) = f{y)- This still allows for a non-trivial dependence of the 
metric on the bulk coordinate. This is the key point in the solution of 0] to the hierarchy 
problem, for localizing gravity in the bulk [|, |^, and for obtaining a non-compact fifth 
dimension with a conventional Newton's force law P|. For an expanding four- dimensional 
universe, however, we must have a ^ n. The metric is determined by solving Einstein's 
equations 

Gab = Rab — -^QabR = f^ Tab- (3) 

Here A, B = 0, 1, 2, 3, 5, and k^ is the five- dimensional Newton's constant, which is related 
to the five-dimensional Planck scale by M^ = k"^. The components of the Einstein tensor 
Gab with the above ansatz are given in (8)-(ll) of 0. 

The five- dimensional stress-energy tensor Tab is the sum of contributions from the bulk 
and from the two branes. The width of the branes is neglected since it is 0(1/M*), which is 
much smaller than the distance scales of cosmological interest. Thus the stress-energy tensor 
is approximated as 

T^^(x, y) = T^^ix, y) + ^f^5{y) + ^^^5{y - 1/2), (4) 

where &o(t) = &(t, 0) and &i/2(t) = &(t, 1/2). Here T^^ is the stress-energy tensor in the 
bulk, and S'^f (5'^jf) is the 4-dimensional stress-energy tensor of the visible (hidden) brane: 

^Bvis = (P' ~-P' ~-P' ~-P' 0) ^^^ ^Bhid = (P*' ~P*i ~P*i ~P*i 0)- ^t this point the composition 
of the energy density on the walls is completely general. 

Using the above energy-momentum tensor BDL derived an equation for the scale factor 
on the visible brane which is independent of the details of the global solution to Einstein's 
equation^: 

^ + (^)' = -|^r,,-^p(p + 3p). (5) 

ao Vao/ oOq 6b 

The time derivatives which appear are with respect to t, where dt = n(r, 0)(ir, and ao(t) = 
a{T^ 0) is the scale factor on the visible brane. There are several important features of this 
equation. First, the energy density and pressure of the second brane do not appear. This 
follows from the local nature of Einstein's equations |T^ . Second, as noted in ^ , the Hubble 



parameter H = do/ao oc p rather than ^/p as in the conventional cosmology. Finally, this 
equation depends only on qq. In the case that T^ is time-independent (as will be the case 
later), this equation completely determines ao(^)- Therefore it is not necessary to determine 
the solutions to Einstein's equations for the whole bulk. 

^This equation can be derived Q] by first calculating the discontinuities in the derivatives of the functions 
a and n at the position of the branes by matching the iJ-functions in Einstein's equations. This information 
about the discontinuity of the derivatives together with the discontinuity and the average value of the 5,5 
component of Einstein's equation results in (pi). 



We assume that no energy is flowing from tlie bulk onto the brane so that T°^ = 0. 
Then the discontinuity in the Gqs = equation [0 gives the usual conservation of energy 
condition: 

p + 3^(p + p) = 0. (6) 

Thus for p = wp, p oc a"^*^'""*'^-' as in standard cosmology. The conservation of energy 
equation (P) and the Friedmann-like equation (j^) are the central equations which are used 
in what follows. 

Now we briefly review the solution presented by Randall and Sundrum @, 0. In their 
scenario the two branes have some tension, Vms, Vud, and the bulk contains a cosmological 
constant A. In terms of our earlier notation, p = —p = Vms, p* = —p* = Vhid, and 
T^ = A(l, 1, 1, 1, 1). A static Lorentz invariant solution is obtained with a{T, y) = n{T, y) = 
e"'^'^\ and b{T,y) = b{y) = Bq = const., the latter having been obtained by a coordinate 
transformation on y. Then the 5,5 component of Einstein's equation gives 

a'' = -^A. (7) 

6 ^ ^ 

Thus A < is required. It is then convenient to introduce 



m — 
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2 - --A. 



Therefore a = ±mbo\y\. (The absolute value is required so that a" is singular at y = 0, 
matching the 5-function sources in Einstein's equations.) The discontinuity equation for 
either a or n then requires that 
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Vvis = T^m = -Vhid- (9) 

We note that the case in which Vms 7^ —Vhid leads to the brane inflating solutions found in 
1^. Since we are interested in non- inflationary solutions, we assume that the tensions of the 
two branes are adjusted to satisfy (P). Since sgniVms) is crucial in determining the expansion 
rate of our wall discussed below, we emphasize that the correlation between the sign of the 
tension of the visible brane and the growth of the scale factor away from the visible brane 
is: 

RSI : niy) = a{y) = e+'-^°l^l ^^ r„. < 0, (10) 

RS2 : n{y) = a{y) = e"™^"!^! ^^ y„. > 0. (11) 

For RS2, gravity is localized about our brane p. This allows for a non-compact fifth di- 
mension that is consistent with the short- distance force experiments. By contrast, case RSI 
with Vms < provides a potential solution to the hierarchy problem [^. The reason for this 
is the following. The metric on the distant brane contains a conformal factor e'^''° in units 
where the conformal factor on our wall is one. This implies that mass scales on our wall and 
the distant wall are then related by niyis = e~"^''°^'^mhid, so that a large hierarchy of mass 
scales is possible if rrihid ~ Mpi ~ M* and mbo ~ 100. 



This non-trivial scale factor significantly modifies the naive relation Mpi ~ M^bo between 
the fundamental and derived Planck scales. The correct relation for m6o ^ 1 is in fact 0] 

1 k'^ 

SttGn = 772~ = '^^"^ = ^\Vvis\- (12) 

Mpi 6 

It is remarkable that this is independent of the size of the extra dimension. For RSI the 
coupling of the KK excitations of the graviton to matter are given by ~ e"^^°^'^/m [Q. Since 
these excitations would appear as resonances in collider experiments, this coupling must be 
0(TeV~^) or smaller, thus implying m ~ Mpi (since e"^'^"''^ ~ 10^^). Therefore a satisfactory 
resolution to the hierarchy problem requires m ~ k~^'^ ~ Mpi. 

However, both of these solutions to Einstein's equations are static and do not describe 
a universe with time-dependent scale factor a. While some work has been done on inflating 
solutions P] , we will attempt to uncover cosmological solutions which reproduce the successes 
of the usual Friedmann equation for a flat universe. This seems difficult at first thought, 
given the earlier statement that in five-dimensional brane models H oc p. However, we will 
see that the presence of large background cosmological constants changes this conclusion, 
and the usual Friedmann equation is reproduced, up to the sign of the source term. 

We begin by perturbing the RS solution by placing an additional energy density on the 
two branes without a compensating change to the bulk cosmological constant. That is, 
we consider p = Kis + Pvis and p = -Kis + Pvis where Ki^ is given by (|) and p^is, Pvis 
are the energy density and pressure measured by an observer living on the visible brane, 
with equation of state Pms = wpms- Key to our results will be that we work in the limit 
Pvis ^ |K>is|- Given that \Vms\ ~ ^p; in these models, this limit is the correct one for 
describing our (post-inflationary) universe. Substituting these expressions for p and p into 
(|) gives, for either RSI or RS2, 

•■ / • \ 2 4 4 

r I I =^ ';r^^vis\'~>Pvis Pvis) TTT^PvisyPvis ~r oPms)- \^'^) 

Oq Vao/ oD 6b 

The 0(A) and 0(/€^V^^^) terms cancel using (||) and (|^). Note that the presence of the back- 
ground energy density allows for H^ oc p^ig as in conventional cosmology. This differs from 
the observations of BDL because here the brane matter is a perturbation to the background 
RS solution. It is also clear that the presence of the prefactor Vyis in ([T3|) implies that for one 
of RSI or RS2 solutions, (^) will have a negative sign relative to the conventional Friedmann 
equation. In fact, we find that the "wrong-signed" Friedmann- like equation corresponds to 
RSI, the solution with Vy^g < 0. To see this, substitute the formula for the Planck mass, 
(ig), into dg): 

RSI: HI—) = i^Pvis- Pvis) - TZ? Pvis (Pvis + ^Pvis)- (14) 

For small densities {pvis ^ Mpi) we can neglect the second term on the RHS, since k^ ~ 
1/Mpi. The first term on the RHS contains a negative sign relative to the Friedmann equation 
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in the conventional cosmology. Obviously this sign is flipped in the RS2 scenario and we get 
the correct Friedmann equations up to small corrections. 

It is interesting to note that for a radiation-dominated (RD) universe this sign problem 
has no effect: since pms = Pvis/'i the first term on the RHS of (|14D identically vanishes, and 
one obtains the same equation for both RSI and RS2. It is then important to see whether 
the conventional RD solution is obtained in this case. The approximate solution for the scale 
factor can be easily found. We take 



ao(t) 



-4 



then 



Pvis{t) = pvtsiU) — 7-T , (15) 

XO-oKh) J 



(16) 




Thus at late times this reduces to the evolution of standard RD cosmology. Note, that 
this solution is different from the ao(t) ~ t^^^ solution presented in (even though (|) 
reduces to exactly the same equation in the case of a RD universe). The reason is simply 
that a non-linear second order differential equation can have more than one solution, and 
the initial condition will determine which is the relevant one. Since ( plGf) reproduces the 
standard RD cosmology, it is plausible that the potential problems with BBN found in |0 
are in fact solved by the existence of the extra solution (|16|). From the solution ([l6|) we can 
see that both RSI and RS2 reproduce the conventional RD solution. Since for RS2 (|T3]) gives 
the conventional Friedmann equation for any type of matter (up to small corrections), we 
conclude that the RS2 solution is viable. Therefore in the following we will only concentrate 
on the RSI solution. 

What is the effect of the wrong sign in the Friedmann equation for the RSI solu- 
tion? Assume that after some time teq the energy density of the universe is dominated 
by a component with an equation of state having w < 1/3. In what follows, the subscript 
"eg" will denote quantities measured at time t^g of matter-radiation equality. Then, e.g., 
Pvis{t) = Peqia-Qit) / cio{teq))'^^^^^^ ■ Ncxt, for w ^ 1/3 and energy densities p^s ^ Mpii the 
second term on the RHS of (IT^) is subdominant to the first, so it is neglected. It is then 



convenient to introduce new variables t = H^qt and x{t) = ao{t)/ao{teq), where H{t) is the 
Hubble parameter. The initial conditions at t = t^q are then x{tf.q) = 1 and x{teq) = 1, where 
the overdot denotes a derivative with respect to i. Then in these units (0) is 

X /x\^_ A 1 

x^\x) ~~2;^^3(^^^' y^^) 

where the dimensionless constant A = 8nGNPeq{^ — 3ti')/(3ifgg) is introduced. Note that 
for w < 1/3 one has A > 0. (One obtains conventional cosmology by the replacement 
A -^ —A.) Also note that if H^ = SttGnp/S were the correct flatness constraint equation 
then A = 1 — 3w. Since it is expected that H^ ~ G^p, then A ~ 0(1). Multiplying the 
above equation hj y = x"^ then gives 

y = _Ay-(3-+l)/2. (18) 
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That is, the expansion of our universe is described by the one-dimensional motion of a 
particle in the classical potential 

V{y) = -^Xy^'-'-^'\ (19) 

For w < 1/3 this describes an attractive potential. From this we understand that the universe 
reaches a maximum size, and then collapses. This is one of the main points of this letter: this 
result is in conflict with our current understanding of cosmology from several viewpoints. 
At the very least, observationally the universe seems to be accelerating [|I^ rather than 



decelerating. (We will examine a sharper conflict with BBN and the age of the universe 
below.) By contrast, in conventional cosmology there is an extra negative sign on the RHS 
of (pISD , so for w < 1/3 the potential is inverted and the universe expands forever (since 
curvature terms have not been included). In the RSI scenario, the period of this oscillation 
is 0(1) in these dimensionless units if A ~ 0(1). In fact, for w = (non-relativistic matter), 
the solution is 

xdx 2 / 1 + A — Ax 



VT+X^^Xx 



1 - A Vl + A - Xx. (20) 



In RSI, the period tu of oscillation (i.e., the age of the universe) for A = 1 is then tu = 20/3. 
So in the original units tu ~ H~^ as expected. By inspection the maximum size of the 
universe is x^ax = 1 + 1/A for w = 0. Finally, note that for the conventional MD cosmology 
(A = —1) the above formula correctly reproduces t ~ x^'^. 



In the RSI scenario, the classical potential (|T^) determines a first-order equation for x, 
in other words, it determines H^ = H^ x^ /x"^. In fact, 

rr2 _ ^^eq SttGn . . 

~ ~2^^ —Pv^s■ (21) 

This contains an arbitrary integration constant E = 1 ^ 2A/(1 — 3w). Note that the above 
equation (^) together with the conservation of energy (|^) implies the Friedmann equation 



(0) for arbitrary values of the integration constant E. In conventional cosmology A 



— 1 -|- "iw so E = 0, and there is an obvious extra minus sign in the second term, reducing 



T|) to the usual flatness equation. In RSI (21) is the analog of the flatness equation. In 



what follows, we assume that p^is > 0, which implies 



3H? 



E > -—^ > . (22) 

eg 



The preceding arguments indicate that in the RSI scenario, after the time of matter- 
radiation equality, the universe collapses on a timescale given by H~g^. The expansion rate 
Heg is obtained as follows. Assuming standard BBN, we know the radiation temperature 
Tbbn ~ MeV, and expansion rate Hbbn during BBN. We use the present-day values of 
the radiation and baryon energy densities, together with conservation of energy and (|16|) 
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describing the RD cosmology, to determine Heq. We use a ~ t^/^ from (|I6|) rather than the 
solution a ~ t^/^, which is known to have difficulties with the He^ abundance [0. Then 



Teq ClBBN ( tsBN] f Heq ^ 



TbBN deq \ teq J \HbBN 

The first equality follows from energy conservation, and the last two equalities from (IT 



(23) 



Next, conservation of energy is used, together with the present-day value p-y/pcrit ~ 10 ^, 
and Tnow ~ 2.7 K, to give Teq ~ 5 eV. Inserting this result, together with Tbbn ~ MeV, 
and Hbbn ~ T^bn/^pi ~ IQ-^i MeV, into (|2|) gives Heg ~ lO'^^ MeV (the standard BBN 
result). Therefore, a standard BBN cosmology implies in the RSI scenario that an MD 
universe collapses on a time scale of tu ~ H~g^ ~ few x 10^ years. It is clear from these 
arguments that the age of the universe has not been used as an input, so this last result may 
be viewed as the RSI standard BBN cosmology prediction for the age of the universe. 

In order to evade the previous arguments RSI requires some non-standard version of 
nucleosynthesis. In particular, in order for the universe to exist for billions of years, the 
Hubble parameter at the start of MD must be orders of magnitude larger than in standard 
BBN:i7f^~103x(fG'^p„,). 

To conclude, we have considered the cosmology of five-dimensional theories with local- 
ized gravity on a three-brane introduced in [^, ^. We have found that the solution with a 
non-compact extra dimension is potentially viable, since it reproduces the conventional cos- 
mological solutions. However, the solution which may solve the hierarchy problem predicts a 
contracting universe with a lifetime of a few thousand years (assuming standard BBN). One 
way to avoid this difficulty might be a non-conventional BBN. But it seems more likely that 
these problems could be avoided by introducing additional matter fields in the bulk (which 
are anyway required to stabilize the radius of the extra dimension). The difficulty with this 
possibility will be maintaining the features which led to the solution of the hierarchy problem 
in the ffist place. 
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